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We discuss the bounds on the mass of Dark Matter (DM) particles, coming from the analysis 
of DM phase-space distribution in dwarf spheroidal galaxies (dSphs) . After reviewing the existing 
approaches, we choose two methods to derive such a bound. The first one depends on the information 
about the current phase space distribution of DM particles only, while the second one uses both 
the initial and final distributions. We discuss the recent data on dSphs as well as astronomical 
uncertainties in relevant parameters. As an application, we present lower bounds on the mass of 
DM particles, coming from various dSphs, using both methods. The model-independent bound 
holds for any type of fermionic DM. Stronger, model-dependent bounds are quoted for several DM 
models (thermal relics, non-resonantly and resonantly produced sterile neutrinos, etc.). The latter 
bounds rely on the assumption that baryonic feedback cannot significantly increase the maximum 
of a distribution function of DM particles. For the scenario in which all the DM is made of sterile 
neutrinos produced via non-resonant mixing with the active neutrinos (NRP) this gives m NRP > 1.7 
keV. Combining these results in their most conservative form with the X-ray bounds of DM decay 
lines, we conclude that the NRP scenario remains allowed in a very narrow parameter window 
only. This conclusion is independent of the results of the Lyman-alpha analysis. The DM model 
|' in which sterile neutrinos are resonantly produced in the presence of lepton asymmetry remains 

viable. Within the minimal neutrino extension of the Standard Model (the t/MSM), both mass and 
the mixing angle of the DM sterile neutrino are bounded from above and below, which suggests the 
possibility for its experimental search. 
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The nature of Dark Matter is one of the most intriguing questions of particle astrophysics. Its resolution would 
Q\ ■ have a profound impact on the development of particle physics beyond the Standard Model. 

Although the possibility of having massive compact halo objects (MACHOs) as a dominant form of DM is still 
qq , under debate (see recent discussion in |l[ and references therein) , it is widely believed that Dark Matter is composed of 
non-baryonic particles. However, the Standard Model of elementary particles does not contain a viable Dark Matter 
particle candidate - a massive, neutral and long-lived particle. Active neutrinos, which are both neutral and stable, 
■ form structures in a top-down fashion [2, La, S 111, [g] , ano - thus cannot produce the observed quantity of early- type 
galaxies [see e.g. 0, H|- Therefore, the DM particle hypothesis implies the extension of the Standard Model (SM). 

The DM particle candidates may have very different masses (for reviews of DM candidates see e.g. [9l. [Tol. TTTL fl^| ) : 
massive gravitons with the mass ~ 10~ 19 eV [13j] . axions with the mass ~ 10 -6 eV [Til, sterile neutrinos having mass 
in the keV range [l5| , sypersymmetric (SUSY) particles (gravitinos [l6[ , neutralinos d7| , axinos with their masses 
ranging from eV to hundreds GeV, supersymmetric Q-balls [Hj], WIMPZILLAs with the mass ~ 10 13 GeV [2(1 l2ll| . 
and many others). Thus, the mass of DM particles becomes an important characteristic which may help to distinguish 
between various DM candidates and, more importantly, may help to differentiate among different models beyond the 
SM. 

It was suggested in ;22j that quite a robust and model-independent lower bound on the mass of DM particles 
can be obtained by considering phase space density evolution of compact astrophysical objects, most notably dwarf 
spheroidal satellites (dSphs) of the Milky Ways. The idea was developed further in a number of works (see e.g. [23|, 

ii MM S3). 

Another way to distinguish between various DM models, and particularly to put a bound on the DM mass, is the 
analysis of the Lyman-a (Ly-a) forest data [2^, [33, [3lt| - 1 This method essentially constrains the possible shape of the 
power spectrum of density fluctuations at comoving scales ~Mpc. Assuming a DM model, (i.e. a particular primordial 
velocity distribution of DM particles), one can obtain a relationship between the DM particle mass in this model and 



1 Absorption feature by neutral hydrogen at A = 1216 A at different redshifts in the spectra of distant quasars. 
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the shape of the power spectrum, probed by Ly-a. 

Although very promising, the Ly-a method is very complicated and indirect. First of all, under the assumption that 
the distribution of the neutral hydrogen traces that of the DM, one can reconstruct the power spectrum of density 
fluctuations at redshifts z ~ 2 — 5 from the statistics of Lyman-a absorption lines. One can then perform a fit of the 
Lyman-a data (often together with the measurements of anisotropy of temperature of cosmic microwave background 
and the data of large-scale structure surveys), to extract the information about cosmological models. This is usually 
done by using the Monte-Carlo Markov chain technique [3^|. At redshifts probed by Ly-a, the evolution of structure 
has already entered the (mildly) non-linear stage. Therefore, to properly relate the measured power spectrum with the 
parameters of a given cosmological model one would have to perform a prohibitively large number of hydrodynamic 
numerical simulations. Therefore, various simplifying approximations have to be realized [30l. l33l. l34l. l35l. l36l. l37l l38l 1391] . 

Apart from these computational difficulties, the physics entering the Ly-a analysis is complicated, and not yet 
fully understood (see e.g. [H EE S3, SI S3 ) . Moreover .the DM particles can significantly influence the background 
physics, further complicating the Ly-a analysis [HI, EH, E?]] . For a recent overview of the Ly-a method see e.g. [48| . 

The systematic uncertainties associated with both computational difficulties and complicated physics of Ly-a sys- 
tems are not fully explored. Therefore, it is very important to have a lower mass bound on DM particles from more 
direct and simple considerations. In this paper we discuss the Tremaine-Gunn and related DM mass bounds based on 
the phase-space density considerations as well as possible ways to strengthen them for several DM models. The ob- 
tained phase-space density bounds are weaker yet comparable with Ly-a bounds and therefore provide an interesting 
alternative. We consider a class of the so-called "generic" DM models, where DM particles are produced thermally and 
decouple while being relativistic, thus having the (relativistic) Fermi-Dirac momentum spectrum. We also consider 
models of non-thermal DM production. In this case the primordial velocity spectrum of DM particles depends on the 
details of the production mechanism. We analyze the case when the velocity spectrum can be approximated by the 
rescaled Fermi-Dirac spectrum, or has two such components (a colder and a warmer one). 

A very important example of such a DM particle is the sterile (right-handed) neutrino. Although known as a DM 
candidate for some 15 years [l5|, recently sterile neutrinos have attracted a lot of attention. It was shown [S§] that if 
one adds three right-handed (sterile) neutrinos to the Standard Model, it is possible to explain simultaneously the data 
on neutrino oscillations (see e.g. [50, EE EH for a review) and the Dark Matter in the Universe, without introducing 
any new physics above electro-weak scale Myy ~ lOOGeV. Moreover, if the masses of two of these particles are between 
~ 100 MeV and electro-weak scale and are almost degenerate, it is also possible [53] to generate the correct baryon 
asymmetry of the Universe (see e.g. [HI, EH)- The third (lightest) sterile neutrino can have mass in keV-MeV range 2 
and be coupled to the rest of the matter weakly enough to provide a viable (cold or warm) DM candidate. 

This theory, explaining the three observed phenomena "beyond the SM" within one consistent framework, is called 
the vMSM [HHl (see also [HI). 

Although weakly coupled, the DM sterile neutrino in the z/MSM can be produced in the correct quanities to account 
for all of the DM. There are several mechanisms of production: non-resonant active-sterile neutrino oscillations (non- 
resonant production mechanism, NRP) [HI l62l l63l [6^. [65|. resonant active-sterile neutrino oscillations in the 
presence of lepton asymm etry (r esonant production mechanism, RP) [H, EH, E3, [HI, decay of the gauge-singlet 
scalar field [6i| (see also [57l . [Tol . l7l|). The Ly-a analysis of the sterile neutrino DM, produced via NRP scenario 
was performed in a number of work s [72l [T3L l74j . These bounds were recently revisited in [48|, using the SDSS 
Ly-a dataset together with WMAP5 [75|. The lower bound on the DM mass was found to be in this case 8 keV (at 
99.7% CL). Ref. (48| also analyzed a more general case CWDM case : a mixture of NRP sterile neutrino with cold 
DM (see also [zij). These results were applied to the RP produced sterile neutrino in [77]]. It was shown that the 
mass as low as 2 keV is compatible with Ly-a data. 

In this paper we will analyze in detail restrictions on the sterile neutrinos, produced via first two production 
mechanisms and briefly comment on the third one in the Discussion. In the case of non-resonant production the 
primordial velocity spectrum is approximately proportional to the Fermi-Dirac distribution [6211781] (the exact spectrum 
was calculated in [m.[65|). 

The paper is organized as follows. In Section |TT] we review DM mass bounds, based on the phase-space density 
arguments. In Section[lTT]we introduce the concept of maximal coarse-graining and propose a conservative modification 
of the original Tremaine-Gunn bound. In Section ITVl we analyze new observational data on recently discovered dSphs 
(see [79, 80] and references therein), and use it to determine the phase-space density of these objects. Special attention 
is paid to determine various systematic uncertainties of measured values. Our results are summarized in SectionfV] We 
conclude with the discussion of the results, analysis of possible uncertainties and outlook for the further improvement 
of the mass bounds in Section IVI1 



2 There are several interesting astrophysical applications of keV sterile neutrinos (see e.g. [45> . 56, 57, 58, 59, gfj and references therein). 
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II. DM MASS LIMITS 



If the DM particles are fermions, there is a very robust bound on their mass. Namely, due to the Pauli exclusion 
principle, there exists the densest "packing" of the fermions in a given region of the phase space. Decreasing the 
mass of DM particles, one increases the number of them in a given gravitationally bound object, containing DM. 
The requirement that the phase-space density of the DM does not exceed that of the degenerate Fermi gas leads 
to the lower mass bound. For example, for a spherically symmetric DM-dominated object with the mass M within 
the region R, one obtains the lower bound m DEG on the DM mass by demanding that the maximal (Fermi) velocity 
of the degenerate fermionic gravitating gas of mass M in the volume |7rii 3 does not exceed the escape velocity 

_ (2G N M\ 1 /' 2 . 
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Here and below g denotes the number of internal degrees of freedom of DM particles, and Gn is the Newton's constant. 
Such a consideration, applied to various DM dominated objects, leads to the mass bound, which we will call m DEG in 
what follows (see Table [TTI below) , 3 

The above considerations assume that the dSphs are purely spherical systems. Analysis of [§3[ shows that ellipticity 
of stars in dSphs vary from 0.22toi2 for Leo IV to 0.80 ± 0.04 for Ursa Major I. Simulated DM halos on the other 
hand tend to have rather moderate ellipticity, cdm ^ 0.32 [H| . 4 According to Appendix [Al the ellipticity of DM 
halos can lower the resulting limit on m DEG by < 10%. 

The limit, obtained in such a way, is very robust, as it is independent of the details of the formation history of 
the system. The only uncertainties associated with it are those of astronomical nature: systematic errors in the 
determination of velocity and density distribution. All these issues will be discussed below (Section IIV1 IIII[) . 

For particular DM models (with the known primordial velocity dispersion) and under certain assumptions about 
the evolution of the system which led to the observed final state, this limit can be strengthened [22|, [H, [24], [2j| Hf| 
HH US HI]- The argument is based on the Liouville's theorem (see e.g. [86|, H^|) and assumes that the collapse of 
the system is disipationless and collisionless. The Liouville theorem states that the phase-space distribution function 
f(t, x, v) does not change in the course of disipationless collisionless dynamics. The consequence of the Liouville 
theorem is that the function f(t, x, v) "moves" in the phase-space, according to the Hamiltonian flow, and therefore 
its maximum (over the phase space) remains unchanged. Therefore, if one could determine the characteristics of a 
phase-space distribution function from astronomically observed quantities (in the first place average density p and 
velocity dispersion a) 5 in dSphs (or any other DM dominated objects), the Liouville theorem would allow to connect 
the measured values with the primordial properties of DM particles. 

One such characteristics of the phase-space distribution is its maximum. Any physical measurement can probe only 
the phase-space distribution, averaged over some phase-space region - a coarse-grained phase-space density (PSD) (as 
opposed to exact or fine-grained PSD). Such a coarse-grained PSD, averaged over phase-space cells AU(x, v) centered 
around points (x, v) in the phase space, is defined via 

f(t,x,v)= / dIL'f(t,x',v') (2) 



voi(An) 

An(i,t) 

(here vol( All) is the volume of the phase-space cell) . From the definition @ it is clear that the maximal (over the 
whole phase space) value of the coarse-grained PSD f max (t) cannot exceed the maximal value of the corresponding 



3 The spatially homogeneous DM distribution is only an approximation. In reality one should consider self-gravitating degenerate fermionic 
gas. It is possible to show that, under some external conditions, the system of weakly interating fermions undergoes a first-orger 
phase transition to a nearly degenerate "fermion star" |8l| . The existence of such objects may also have insteresting astrophysical 
applications |82( . 

4 Therefore it is hard to explain the ellipticity of stars in the most elongated dSphs, see the discussion in |83| . 

5 The quantity directly measured in astrophysical observations is the projection of stars' velocities v(R) along the line of sight. We will 
denote such a projection by V(R) to distinguish it from the absolute value of the 3D velocity v(r). The ID velocity dispersion is 
denned as a(R) = <y 2 (fi)> 1 /2 an d is in principle the function of the projected radius R. However in the DM dominated objects, for R 
greater than certain characteristic scale rotation curve flattens, a const. It is this constant which is usually referred to as "velocity 
dispersion" . 
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fine-grained PSD. On the other hand, as a consequence of the Liouville theorem, the maximum of the fine-grained 
PSD f m ax does not change in time. Thus, one arrives to the following inequality 

fmax(fy ^ fmax ■ (3) 

The inequality §5§ allows to relate the properties of DM at present time t with its primordial properties, encoded in 
fmax- For example, if one assumes that initially DM particles possess relativistic Fermi-Dirac distribution function 
with some temperature T FD (relativistically decoupled thermal relics): 

/fd(p) = (2^)3 + 1 (4) 

and recovers from astronomical measurements that in the final state the coarse-grained PSD of the system is described 
by the isothermal sphere (see e.g. (89l |) with a core radius r c and a ID velocity dispersion er, whose maximum is given 
by 

9cr 2 

fiso,max = ~. "^J TZ 2 (^) 

the comparison of the maximum of the coarse-grained PSD ([S]) with its primordial (fine-grained) value leads to the 
so-called Tremaine-Gunn mass bound [22I ]: 

m „> mTO , where <, s „ **f (6) 

(2tt ) b / 2 gG N a r 2 

For the case of initial distribution |[4]) this bound is stronger than the one, based on the Pauli exclusion principle, by 
a factor 2 1 / 4 [22^ . For different primordial DM distributions this difference can be significant (as we will demonstrate 
later). We would like to stress, though, that these stronger bounds make assumptions about the evolution of phase- 
space density, while the one, based on the Pauli exclusion principle does not assume anything about either primordial 
velocity distribution of the particles, or the formation history of the observed object and simply compares measured 
phase-space density with the maximally allowed for fermions. 

Another characteristics of the phase-space distribution function is the "average phase-space density" 

introduced in [26|, [27j. The value of Qf (average PSD today) is simply defined in terms of the observed quantities p 
and (v 2 ) = 3cr 2 and therefore serves as a convenient estimator of the PSD for any DM dominated object. One can 
calculate primordial Qi for an arbitrary homogeneous distribution function f(p) 
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and compare it with its value today Qf. It was claimed in [261 |27j that Q cannot increase during the evolution of 
DM: 

Qi > Qf- (9) 

Applying this inequality to the dSphs, one obtains several times stronger mass bound, than that of [22l |. 

To illustrate the origin of the inequality ([5]), authors of [26l.[27| noticed that in the case of the uniform monoatomic 
ideal gas, Q is related to the usual thermodynamic entropy per particle (see Appendix IB 1[) and the inequality for Q 
becomes a consequence of the second law of thermodynamics. Indeed, in this case one can see that 

M._ l0g («££5!) +l06C[/1 , (10) 



where in the right hand side of (TIT))) functional C[f] does not depend on the average density and velocity of the DM 
particles. 
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However, because of the long-range interaction of DM particles, the notion of Boltzmann entropy is well-defined 
only for the primordial DM distribution and not for the final state of DM evolution (see e.g. the discussion in [90(). 
Moreover, we will show below that in general the increase of entropy does not imply the decrease of Q. Indeed, the 
values of C\f \ are different for different types of phase-space distributions / and therefore they can change with time 
if the shape of the (coarse-grained) distribution changes. Namely, even if initial (i) and final (/) states both satisfy 

relation (flU)) between the entropy and Q [Si / = logQ t — log Qi '{ h ) from the second law of thermodynamics 

S f > Si (11) 

it only follows that 

Q*>Qf^- (12) 

Therefore, in general, the inequality (J9j) does not follow from entropic considerations. 

Moreover, the simple relation (|10|) between the entropy and Q does not hold for the distributions we are interested 
in. For example, for the Fermi-Dirac distribution Q one has: 



S_ 

N Qh 



const^ (13) 



(see Appendix IB 31 for details). The relation becomes even more complicated, if one considers DM candidates (e.g., 
sterile neutrinos, gravitinos), which are produced out of thermal equilibrium. In general, when the primordial dis- 
tribution function depends on several parameters, both Q and entropy are expressed through these parameters in 
a non-trivial way and the simple relation (|10[) does not hold. For example, this is the case when DM is produces 
in two stages and the DM distribution shape has two components: colder and warmer one. Physically interesting 
examples include: production of sterile neutrino in t he p resence of lepton asymmetry [6(1 l67l . 1(381 ]: production of 
gravitino thermally at high temperatures (see e.g. [91, 92]) accompanied by non-thermal production via late decays 
of next-to-lightest supersymmetric particles (see e.g. (93|| ) . 

Keeping in mind the above considerations, one might be tempted to use the entropy of the system as an estimator 
of PSD and utilize the entropy increase (TlTj) instead of the inequality on Q to put a lower bound on the DM 
mass. However, unlike Q, which by definition is expressed solely in terms of measured quantities p and a, the 
inequality (jlip requires the knowledge of the phase-space distribution function in the final state (e.g. to determine 
the Cf in the right-hand side of Eq. (|10p or, more generally to express the entropy of the final state in terms of the 
observed quantities). This information cannot be simply deduced from astronomical observations. One possible way 
to formulate a conservative, robust inequality would be to find the maximal possible entropy for a given system with 
measured macroscopic parameters. However, it was shown in [86l [891 . l94l [95| that such a maximum does not exist. 
Namely, for a gravitating system which usually consists of a compact core and a widely dispersed halo of finite mass, 
the total Boltzmann entropy of the system goes to infinity when the halo becomes infinite. Physically, the measured 
density and velocity dispersion characterize the inner part of the object. The astronomical observations do not usually 
probe the outskirts of gravitating systems (such as dSphs) and phase-space distributions (such as (0) do not describe 
them properly. On the other hand, to compare with the homogenous initial system having a primordial velocity 
spectrum, we need to know an entropy of the whole system. The large (and unknown!) fraction of this entropy can 
be related to the outskirts. The entropy of the gravitating system depends on the precise state of the halo. 

As a result, it is not possible to construct a simple and robust limit, using entropy considerations. 



III. MAXIMAL COARSE-GRAINING 



In view of the above arguments, to derive a conservative mass bound, in this work we will follow the original 
approach of Tremaine and Gunn [221 ] with some modification. 

An important advantage of this approach is that the maximum of the phase space density is likely to be located 
in the inner, dense part of an object. Therefore, under this reasonable assumption, the results do not depend on the 
DM distribution in the outskirts (see the discussion above). 

As discussed already, the coarse-grained phase-space distribution in the final state cannot be measured directly, 
and one has to make assumptions to deduce its maximum. A conservative way to minimize this uncertainty is to 
use the "maximally coarse-grained distribution" . It is based on a simple fact that the mean value of a function, 
averaged over an arbitrary region cannot exceed its maximal value. Therefore, the average value of coarse-grained 
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phase space density in a large phase-space volume can be taken as a conservative estimate of the F max , independent 
on assumptions about the actual form of phase-space distribution. 

To this end we consider an (approximately spherically symmetric) gravitating system (having in mind a dwarf 
spheroidal galaxy), that has the mass M(R) confined within the radius R. The phase-space volume, occupied by the 
DM particles, forming such a system can be approximated by 

IIoo = (|tt) R 3 vl , (14) 

where we have introduced escape velocity v 2 ^. The "coarsest" PSD is such that the averaging ^ goes over the whole 
phase-space volume: All = IIoo: 

p_ M 9 M 3p 

As an estimate for R we take half-light radius ru (he. the radius where surface brightness profile falls to 1/2 of its 
maximal value). Neglecting possible influence of ellipticit y o f stellar orbits (c.f. Appendix |A|) . assuming constant 
DM density within and isothermal distribution of stars [96| , we obtain the following estimate on the average DM 
density within r^: 

_ 3 log 2 a 2 

P=^ n — 2> ( 16 ) 

2?r Gjvr£ 

Assuming isotropic velocity distributions, 6 the escape velocity Voo of the DM particles is related to the velocity 
dispersion a via ~ y/6cr. In such a way we obtain the averaged PSD F: 

p _ M _ P M 31og2 R , 125 M Amy 3 / km/sec \ flpcV 
Tl x 8ttV6ct 3 ~ 16V6ir 2 G N ar 2 h pc 3 V sec/ \ a ) \ r h ) ' 

which coincides with its maximal value (being flat). 

As a consequence of Eq. ([3]), this "coarse-grained" PSD F is smaller than the f m ax ~ the maximum value of 
fine-grained PSD, equal to its primordial value: 

F < f ma x ■ (18) 

Eq. (|18p relates the observed properties of the DM-dominated systems (l.h.s.) with the microscopic quantity on the 
r.h.s. of inequality, which depends on the production mechanism of the DM. 

In this paper we are mostly interested in two types of primordial momentum distribution. One is the relativistic 
Fcrmi-Dirac @i with its f m ax being equal to 

g 

fmax ' FD = 2(27r!T) 3 (19) 

(we fix the overall normalization of the phase-space distribution function by the relation M — J d 3 x d 3 v f(t,x,v), 
where M is the total mass of the system). Another one is an (approximate) form of the momentum distribution for 
sterile neutrinos, produced via non-resonant oscillations with the active ones [HI, |62[. F° r t ne latter case we consider 
the velocity dispersion to be 7 

/-b) = ^vfri ■ (2°) 



This assumption seems to be correct for the DM particles, since numerical simulations of DM structures of different scales show that 

CT 2 +(T 2 

the velocity anisotropy j3(r) = 1 ^ * tends to be zero towards the central region [97l, |98|, l99l llOCt llOll . It is not clear whether (3 

equals to zero for stars in dSphs. The assumption of isotropy of stellar velocities leads to the cored density profile s |l02j, | l03 |. therefore 
our estimate for p tends to be robust. This is confirmed by comparison of the estimate H16H with those, based on 104, 105, 106), where 
DM density profiles were obtained under the assumptions of different anisotropic distributions of stars in dSphs. 

In reality the momentum distribution in the case of non-resonant production does not have thermal shape. The exact shape, taking 
into account contributions from primeval plasma at temperatures around QCD transition, can be computed only numerically [64L l65|. 
The difference between the exact distribution and (I20H does not exceed 20%, which does not affect the mass bounds. 
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The normalization constant \ lS proportional to the mixing strength between active and sterile neutrinos and T v is 
the temperature of neutrino background T v (z) = (l + z)T U0 , related to the temperature of the CMB background today 
via T„ = (4/ll) 1 / 3 T CMB .o. For the maximal value of distribution (f20|) we find 



gxm 
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f » A. " NRP fOT \ 

From the definition ([20]) one can relate the normalization factor g\ to the DM abundance (see e.g. fF^ ) 



_ n ,2 m NRP [eV] 
w DM = S2 DM ft = ffX 94 eV • (22) 



Therefore we can rewrite maximal value of the primordial phase-space density (|22l) as 



2(27rfr) 3 eV 3 



_ ^ DM '"NRP frjo'l 
Jmax, NRP — o/0_t;\3 --ir3 ' \ ^/ 



Notice, that unlike the Fermi-Dirac case, for the NRP scenario f m ax behaves as the third power of particle's mass. 

In the presence of lepton asymmetry in primeval plasma the resonant production of sterile neutrinos becomes 
possible [601 . A possible lepton asymmetry, generated in the framework of the z^MSM and spectra of sterile neutrino 
DM were recently computed in [63, [fil]- Qualitatively, these spectra contain a "cold" (resonant) component and 
a "warm" one, produced through non-resonant oscillations, analogously to the NRP scenario of [la]. The spectra 
as a whole become colder than in the NRP case (see e.g. Fig. 6 in (68|). The maxima of primordial phase-space 
distributions for these spectra are higher (sometimes significantly) than for spectra, produced in the NRP scenario 
(c.f. Fig. 5 in [68]). Therefore, in general mass bound for such a DM is expected to be weaker than that of the NRP 
scenario. The exact form of these spectra can be computed only numerically. We used a number of spectra 8 to check 
those which satisfy the bound (|18p or TG bound (see Section W} . 

Let us compare expression (|17[) with the original Trcmainc-Gunn bound (maximum of the right hand side of Eq. ©): 



^tg = 1= • (24) 



The values of F is smaller than F TG by 



F log 2y/ir ( rv \ 2 _ _ _ ( r r x " 



_ , , -0.118 — , (25) 
F TQ 6V3 WJ Vh) ' 

where r c and are the core radius of isothermal profile and the half-light radius, correspondingly. When comparing 
F and F TG below, we take ~ r c . Essentially, the difference between F and F TG is due to the different assumed 
velocity distributions. While the Maxwell distribution was assumed in [22| (c.f. Eq. (J5J), we assume constant velocity 
profile from escape velocity Voo down to v = (as shown on the Fig. [T|). The numerical factor in (f2"5]) is the ratio of 
areas under two velocity curves of Fig. [T] Translated into the mass bound, relation (|25[) means that for DM particles 
with distribution |4} one would obtain roughly 40% stronger mass bound by using the original Tremaine-Gunn bound, 
rather than F (and ~ 60% stronger mass bound for the case of the distribution (120|)). 

Let us compare our new bound with the one, based of [26l[27|. Following the definition ([7]), we express the measured 
value Q / for a dSph through the observed quantities 

Q - -^37, - 14.83^% ( *E) 3 f^£l) fill) 2 1 , (26) 
if^a^Y'^ pc^ \ sec / \ a J \ ru J 

where r\ is the scaling factor which accounts for the fact that the dark matter particles do not necessarily have the 
same velocity dispersion as the stars, is the half-light radius, a is the measured one-dimensional velocity dispersion 
of the stars and p is defined in (TT6"]) . It was estimated in [26[ that r/ s» 1 . In Eq. (f2"6")) we used the same value of p as 
in Eq. p?|) . For for the same dSph, Qf is bigger than F (given by expression (|T7|) ) by a factor 87rV2/3 w 11.85 ... . 



We are grateful to M. Laine and M. Shaposhnikov for providing these spectra to us. 
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On the other hand, for any initial momentum distribution f(p) we should compare Qi, given by Eq. p]). with the 
fmax ■ For both types of distribution ((4]) and (|20|) the ratio of initial Qi / P max is given by 

Qi - - 0.973... (27) 



/W a 57T5C 3/2 (5) 

As a result, a bound, based on the decrease of the average PSD Q is stronger than F bound from the same object by 
a factor: 

!™9±n 12.176... (28) 

(where again we put ?y = 1). This leads to ~ 1.87 times stronger bound on the m FD and w 2.3 times stronger bound 
for m NnP . 



IV. ANALYSIS OF MEASURED VALUES 

Recently, a number of very faint, very dense dSphs were detected [zl, M IM EM EH EH E3 • To calculate the 
mass limits, we used the data from two recent papers: [zH H(|. First of all, we should notice that although both of 
these pap ers provide the estimate of Q for each object, they use different prescriptions for computing this value. 

In [80| | the quantity Q is estimated inside the half-light radius r/i, using one-dimensional velocity dispersion a of 
stars: 

Qgu = 4 = a 2 • ( 29 ) 
Compared to our definition (TTT)) F — Qgu ~ O.OAbQcu- Following [l!2i | the authors of [zH define central density 

p = imaW/rl (30) 



where n ~ 1 is a numerical parameter, characterizing plausible density profiles (for details see [79l. |112| ). They used 
Po to define the quantity: 

Qsg = ^ , (31) 

As a result for the same object Qsg is by a factor of 14.60 greater than Qgu- 

Using the available information about dSph galaxies (refs. (79l.[80T| and refs. therein), we calculate F, trying also to 
estimate the errors. Several factors contribute to the errors of a and r^. 
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Galaxy 


rn, Plummer 


rh, exponential 


Coma Berenices 


5.0' 


5.9' 


Canes Venatici II 


3.0' 


3.3' 


Leo IV 


3.3' 


3.4' 


Hercules 


8.0' 


8.4' 



TABLE I: Uncertainties of determination of half-light radius for several dSphs. 







rh 


a 


F 




m FD 




?TInrp,tg 


dSph 


References 


pc 


km/s 


Mq pc -3 ( km/ sec) -3 


keV 


keV 


keV 


keV 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 



dSphsfrom [j(| 



Sextans 


[so, m] 


630±170 


6.6±2.3 


4.78t 8 ,^ 


•10" 


-7 





147 +0.044 

-0.026 





17/1 +0.053 
1 '^-0.031 


U.4O4_ 104 


(.+0.302 
U. 1 10_o.l63 


Fornax 


[so, m] 


400±103 


10.5±2.7 


7.45^ 


■10 


-7 





lfi4 +0.041 
1D4 -0.026 





1Q r+0.049 
lyJ -0.031 


n CTO7+0.183 
U.O^ I -0.108 


0.830lg; 2 f 


Leo I 


[so, m] 


330±106 


8.8±2.4 


i 01+2.59 
1 ' ol -0.72 


■10" 


-6 





18q +0.059 
lo3 -0.034 





99/1 +0.070 
zzy: -0.041 


635+ - 279 
u.ooo_ 14g 


1 00+ 044 

l.UU_ 2 3 


UrsaMinor 


[so, m] 


300±74 


9.3±2.8 




•10" 


-6 





iq 5 +o.05i 

lyJ -0.031 





909+0. 060 
ZJZ -0.037 


0.665l°;?« 


l.UO_ 22 


Carina 


[so, m] 


290±72 


6.8±1.6 


9 1 q+ 2 - 87 


•10" 


-6 





9 i r+0.050 
zld -0.032 





9 rr+0.060 
ZOO -0.039 


0.755l^« 


1 iq+0-38 

l-iy_0.23 


Draco 


[80,83] 


221±16 


9.5±1.6 


2 70+ 107 

z - ' u -0.69 


■10" 


-6 





99fi +0.020 
zzu -0.016 





9 fiQ +0.023 
ZU3 -0.019 


0.809l°;g?» 


1 9 7 +0.15 
L - /L ' -0.12 


Bootes 


[m, 113, 124] 


246±28 






•10" 


-6 





oqc+0.029 
ZOU -0.027 





900+0. 035 
ZOU -0.033 


855+ - 143 
u.ooo_ 130 


1 qc+0.23 
l.OO_ 20 


Sculptor 


[so, m] 


160±40 


10.1±0.3 


4 qq+ 3 -62 


•10" 


-6 





2fi4 +0.038 
zo ^-0.031 





o 14 +0.046 
Jly: -0.037 


u.»yo_ 154 


1 rc+0.312 
l.OU_ 243 


Leo II 


[so, m] 


185±48 


6.8±0.7 


c; qo+5.55 
2.30 


■10" 


-6 





2fiq +0.052 
zua -0.O35 





Q1O+0 062 
J19 -0.042 


1 02+ ' 27 


1 fin+o- 43 

l.DU_ 27 



dSphsfrom [zj 



Canes Venatici I 


[79. 83] 


564±36 


7.6±2.2 


r 1 7 +3.46 
d - 1 ' -1.63 


•10" 


7 


1 ™+ - 020 

U.lOU_ 013 


n 1 78+0.024 

u.ii , 8_ . i 6 


467+0- 087 

_ .055 


r, 7 ot;+0.137 
U. IOJ_o.o87 


Ursa Major I 


[79,83] 


318t§S 


7.6±2.4 


1 fi o+1.46 


•10" 


6 


iqq+0-035 

U-^»»_0 .026 


n 907+0.041 

U.zo 1 _ .031 


0.684«;^ 


1 08+ 26 


Hercules 


[79,83] 


330iS 


5.1±2.4 


9 9 r+3.74 
A.ZO_ 1 2 g 


■10" 


6 


n 21 r+ 060 

U./lO_ 041 


n 9C7+O.O71 
U.ZO^_ .049 


0-762l°; 2 r 7 


1 9n+ - 46 

l.^U_ 29 


Leo T 


[79,83] 


178±39 


7.5±2.7 


■J-^O-2.66 


•10" 


6 


2 67+ 071 

U.ZOf _ 043 


0.318±g;gg* 


l-0llo.ll 


1 5q+o.59 

±.o»_q 33 


Ursa Major IF 


[79,83] 


140±25 


6.7±2.6 


Q irq+13.55 
a ' do -4 .59 


■10 


-6 


31 0+ 077 

u - Jlu -0 .047 


0.369lH 9 56 


1 90+0.42 

± - ZJ -0.24 


1 q4+0-67 


Leo IV 


[79,83] 


116 +_lt 


3.3±2.8 


9 09+34.39 


■10" 


-5 


4nfi +0 - 368 


0.483l°; 4 f 9 


1 77+2.41 
L -l '-0.55 


9 7Q+3.80 
z - ' y -0.87 


Coma Berenices 


[79,83] 


77 ± 10 


4.6±2.3 


4.59t 7 ;g 


•10" 


5 


U.4O»_ 069 


U.O £ lD_o.o82 


9 ns+0- 80 


q 9Q+1-25 
°- zo -0.64 


Canes Venatici II 


[79,83] 


7£\t 


4.6±2.4 


4 q 7 +8.92 

' -2.66 


•10" 


5 


0.468t°; l^ 


fl ccy+0.163 
U.OO 1 -0.097 


9 14 +0.87 
z - i4 -0.48 


Q qR+ 138 



TABLE II: Parameters for dSphs from 79, 80] (columns 1-5) and derived lower mass limits for various types of DM (columns 6- 
9). m DEG refers to the limit from Pauli exclusion principle (fTJ, m FD is the limit for particles with the momentum distribution Q, 
m N Rp and m N Rp,TG - for d istribution (|20p . All results are quoted for g = 2 internal degrees of freedom. Results for NRP scenario 
are for lo dm = 0.105 [l25| . 

"There is an extensive evidence that Ursa Major II is a tidally disrupted dSph [see e.g. |79|| - Therefore, the results for it are provided for 
illustrative purposes only. 

First of all, as a is the dispersion of measured velocities, it has the statistical error (which can be quite large for 
the ultra-faint dSphs where the number of stars can be rather small (~ 10 — 100, c.f. [7^, Table 3]). However, the 
systematic error is much larger. The authors of [79| found the systematic error on their determination of velocity 
dispersion to be 2.2 km/ sec. We add this error in quadratures to the statistical errors, found in [73, Table 3]. The 
results are shown in the column number 4 in the Table [TTJ 

The half-light radius is a derived quantity and there are several contributions to its errors. First of all, the 
surface brightness profile is measured in angular units and their conversion to parsecs requires the kn owledge of the 
distance towards the object. These di s tances ar e gen erally known with uncertainties of about 10% (see jl07l . lll0l . lll3[ 
[Till . [TT1 [HE LUl, Iml . fill . LUl EM LH1 EH UM)- Another uncertainty comes from the method of determination 
of rh- The surface brightness profile gets fit to various models to determine this quantity. For several dSphs: Coma 
Berenices, Canes Venatici II, Hercules and L eo IV authors used two different profiles (Plummer and exponential) for 
evaluating the annular half-light radius [107] . Their results are present in the Table [H We use these results to estimate 
the systematic error on r/j to be 20% and use it for all the dSph, where is quoted without errors. The results of 
determination r^ are shown in the 3rd column of the Table [TTJ The obtained values of F with corresponding errors 
are presented in the Table [TTJ column 5. We determined the errors on F by pushing the uncertainties in both a and 
r^ so that the values of F is minimized (maximized). 
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V. RESULTS 

Our main results are compiled into the Table HT1 (columns 6-9). The column 6 of Table HT1 contains the bound 
on m DEG (given by Eq. JTJ)) based on the Pauli exclusion principle. It is independent of the details of the evolution 
of the system, is not affected by the presence of baryons (see below) and holds for any fermionic DM. The column 
7 contains the mass bounds for the relativistically decoupled DM particles (primordial distribution (|4])), obtained 
by combining Eqs. (fTT)) - (fll?|) . Combining Eqs. (fTT|) , (|18[) and (f2"B")) one obtains the result for the case of DM with 
primordial velocity distribution (|20|. quoted in the column 8. Both bounds in columns 7 and 8 conservatively 
assume maximally coarse-grained distribution function (see Section llllj) . In instead of the maximal coarse-graining, 
one assumes the isothermal distribution in the final state (c.f. Fig. [l}, one arrives to the original Tremaine-Gunn 
bound, shown in the 9th column. It is obtained by comparing the expressions (f2Tj) with (|24|) . 9 We denote the 
corresponding mass bound by m NR p !TG . 

We quote all the mass bounds with the corresponding uncertainties, coming from those of in determination of a and 
Th (see Section llVI) . However, for any given object there can be unique reasons, violating the standard assumptions 
and therefore increasing the uncertainties. Therefore, although the strongest bounds in Table ITT1 come from the Canes 
Venatici II (CVnII) dSph, we decided to take a value which independently follows from several objects as a single 
number, characterizing our results (for a given type of DM). To this end we choose the value, obtained for Leo IV. 10 
Thus, the mass bounds, quoted below are excluded from three dSphs: Leo IV, CVnII and Coma Berenices (Com) 11 
To summarize, we obtain the following lower bounds 

m DEG > 0.41 keV , (32) 

m FD > 0.48 keV , (33) 

m NRP > 1.77 keV , (34) 

and 

m NRP . TG > 2.79 keV . (35) 

We can compare lower bounds (|34|) (|35p wit h th e upper ones, comin g from ast r ophysical (X-ray) c onstraints on the 
possible flux from sterile neutrino DM decay [TH, [l27| [TH [l2l HH HM [TH [lM [TH, (TH, [l36|. Taking central 
value (|34p and comparing it with the X-ray constraints, one sees that there exists a narrow window of parameters 
for which 100% of DM can be made from the NRP sterile neutrino (c.f. Fig|2|). Less conservative bound (|35|) . based 
on [22|] (marked by the dark orange double-dotted vertical line on the Fig. [5]) almost completely closes this window. 
Notice, that these bounds are comparable with the lower mass limit m NRP > 5.6 keV, coming from the Ly-a forest 
analysis of (74J. 

We also performed the analysis for sterile neutrinos, produced in the presence of lepton asymmetry (resonant 
production mechanism, RP) [6y, [67], [68| . This mechanism is more efficient than the NRP scenario and allows us to 
achieve the required DM abundance for weaker mixings (c.f. Fig. 4 in 68]). This lifts the upper bound on the DM 
particle mass in RP scenario to ~ 50 keV. To estimate the lower mass bound at this scenario, we have analyzed a 
number of available spectra (mass range 1 — 20 keV, asymmetries (2 — 700) x 10~ 6 (see [67|, |68[ for the definition of 
asymmetry). The result are collected on the Fig. [3l One can see that based on F, the M RP = 1 keV is ruled out for 
lepton asymmetries L > 10 2 and higher masses M RP > 2 keV are allowed for all available asymmetries. Based on the 
original Tremaine-Gunn bound, Af RP = 2 keV is also ruled out for sufficiently high (L > 10 2 ) lepton asymmetries. 
Thus, resonantly produced sterile neutrinos remain a viable DM candidate (see Fig. 3]). 

Finally, we would like to notice that our bounds (|3"3"|) - ([3"5]) are valid under the assumption that the influence of the 
baryons does not result in the increase of the PSD in the course of structure formation. If this assumption does not 
hold, only the bound (|32[) remains intact. 



The value of r c is not currently known for several new, faint dSphs, from which we obtain the best limits on DM mass. Therefore, to 
calculate the Tremaine-Gunn limit in Table|IT] we use the conservative estimate r c ?s r h (see comment after Eq.((2SJ). 
Notice, that the numbers for Leo IV essentially coincide with the mass limits from CVnII and Com if all uncertainties in these dSphs 
are pushed to minimize the mass bound. 

It is possible that Coma Berenices is undergoing tidal disruption (like another ultra-faint dSph, Ursa Major II (UMall), closely resembling 
Com) [79l | . However, unlike UMall (or the best known example of tidally disrupted dSph, Sagittarius), there are no known tidal streams 
near the position of Coma Berenices and the evidence in favor of tidal disruption are quite moderate [c.f. discussion inl79l. §3.6]. 
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FIG. 2: Restrictions on parame ters of sterile neutrino (mass and mixing sin 2 (28) between sterile and active neutrinos) from 
X-rays f [1281 . [l33l . Il34l . Il35l . Il36l |) and phase-space density considerations (this work). Our analysis excludes the region to the 
left of the vertical line (|34|) (purple shaded region). Two dashed-dotted vertical lines mark the systematic uncertainties of this 
bound. The dotted line on the left marks the bound (J32J based on the Pauli exclusion principle. The double-dotted dark 
orange line marks the bound (|35[l . The black dashed-dotted line is the NRP production curve (i.e. pairs of m mF and 9 that 
lead to the correct DM abundance) 65]. The gray region marked "NRP production" accounts for possible uncertainties in the 
abundance computations (see [63.1651] for details). 



In this work we chose not to use the bound, based on the "average PSD" Q [2|| [27| (see discussion in the Section HI)). 
However, as this bound is widely used in the literature, we quote analogs of lower limits (j33[) and (134| based on 
inequality ^ (which we denote to fd . hd and to NRPjH d correspondingly): 

«i FDjHD = 0.9 keV , 
^nrp,hd = 4.0 keV . 

For details see Appendix [C] 



VI. DISCUSSION 



In this paper we suggested that a conservative way to put the bound on a DM particle mass may be based on the 
requirements that the maximum of the observed coarse-grained phase space density should not exceed the maximum 
of the initial distribution function of the DM particles. The maximum of the coarse-graned distribution function in 
the final state may be conservatively estimated from the observed quantities. This bound relies on the assumption 
that the maximum of the distribution function was not significantly increased by the interaction with baryons. 

Although DM consists of the non-interacting particles, the remaining part of the galaxy - the baryons - interact with 
one another and dissipate their energy, finally concentrating towards the center. The baryons , which ar e condensed in 
the center, influence the shape of DM halo gravitationally, increasing the central DM density [I37lll38j . The opposite 
effect is the energy feedback from SNae, galactic winds and reionization, which creates the strong outflow, significantly 
decreasing the mass of the gas and thereby affecting the DM halo shape. Such a fee dback is thought to be re sponsible 
to the formation of dwarf spheroidals from gas-rich dwarf spiral/irregular gala xies [139L Il40l Il4ll Il42l . Il43| . Clearly 
both gas condensation and feedback strongly influence the central PSD of DM [1441 ] . and in principle can lead to the 
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FIG. 3: Restrictions on resonantly produced sterile neutrinos. Primordial fmax is computed numerically based on the spectra 
fro m [67|. loSll . Different colors parametrize different lepton asymmetries for a given mass (see definition of lepton asymmetry 
in [63, |68|]). Grey shaded region is bounded by the maximal and minimal values of F for Leo IV (from Table ITT1 column 5). 
Horizontal dotted lines represent central value for F (lower) and F TG (upper) for Leo IV. The DM spectrum is ruled out if the 
point falls into the shaded region (below the dotted line). 
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FIG. 4: Allowed window of parameters for sterile neutrinos produced via resonant oscillations (white unshaded strip between 
two black lines). Two bounding black lines are obtained for non-resonant (upper line, lepton asymmetry = 0) and resonant 
production with the maximal lepton asymmetry, attain able in the vMSM [67], EH (lower line). The colored regions in the upper 
right corner represent X-ray bounds |128| . Il33l . Il35l [l3q |. Region below IkeV is ruled out from the PSD arguments (this work). 
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violation of the inequality ([3]). Numerical studies o f gal axy mergers show that baryons can lead to the increase of 
the phase-space density during the merger (see e.g. |145| ). However, the method used in this work - coarse-graining 
of the PSD over a large phase-space region - reduces the influence of baryons. Indeed, we take the spatial averaging 
over the radius R ~ r^, which includes external part of the system, where the amount of baryons is small. Additional 
studies are necessary to estimate effects of baryons and make our bounds more robust. We plan to address these 
issues elsewhere. 

We would also like to stress that the initial velocities of DM particles in our approach are thermal velocities and 
they should not be confused with the so-called Zeldovich velocities |2j . Numerical simulations of galaxy formation do 
not start at the time, when the DM phase-space distribution is spatially uniform (redshifts z > 10 3 ). Instead, the 
initial (linear) stage of the structure formation is computed analytically in the framework of the so-called Zeldovich 
approximation (2j. This approximation is commonly use d to set up initial conditions for the numerical simulations 
of non- linear stage of structure formation [l46l Il47l 1 14=81 ] . which start at redshifts z ~ 10. The peculiar {Zeldovich) 
velocities acquired by DM particles at this stage due to structure formation and included into the initial conditions 
are normally a ~ 10 km/ sec. Apart from Zeldovich velocities, DM particles also possess thermal velocities, which 
are discussed in this paper. For cold enough Dark Matter these thermal velocities are much smaller than Zeldovich 
ones and, thus, are o ften neglected and not includ ed into initial conditions. Therefore, the numerical studies of PSD 
evolution 12 (see e.g. [l49l 1 150L 1 15 ll Il52l . Il53l Il54j ) essentially investigate the change of PSD from Zeldov ich t o final 
stage. It was found in some of these works that the PSD changes by 10 2 — 10 3 in the process of collapse [l50t ] . This 
change of PSD can be understood as being simply an evolution from initial Zeldovich velocities <Zj ~ 10 km/ sec to 
the final (virial) ones 07 ~ 10 2 km/ sec (with Qi/Qf ~ (07/0^ ) 3 ~ 10 3 ). 

Because initial thermal velocities may be much smaller than Zeldovich ones, initial PSD may differ from the final 
(observed) PSD not by 2-3, but by many orders of magnitude. This fact does not contradict to the results of 
simulations, desc ribe d in e.g. [l50l ] and, therefore, cannot be used to obtain an upper bound on the mass of DM 
particles (c.f [H, EH). 

This work was mostly concentrated on restrictions on the mass of the sterile neutrino DM, produced in through the 
non-resonant oscillations with active neutrino (NRP scenario). We see that our results (Section [Vj) strongly disfavor 
such sterile neutrinos as the single DM component. This conclusion is not based on the Ly-ct method and therefore 
is not subject to its uncertainties (discussed in the Introduction). However, several uncertainties can affect this 
conclusion, the major being baryonic feedback. To make this result really robust, apart from further modeling of the 
baryonic influence, one needs to strengthen the tension between upper and lower mass bounds discussed in this paper. 
This is plausible and may be done either by improving the X-ray bounds with new observations or by strengthening 
the PSD consideration, which is in the first place related to better measurements of kinematics of dSphs. 

In the presence of lepton asymmetry, the resonant production (RP) of sterile neutrino DM takes place [66| . This 
mechanism is more efficient [661 . 1 6 71 |68[ than the NRP scenario and allows to achieve required DM abundance for 
weaker mixings (c.f. Fig. 4 in [68]). This lifts the upper bound on the DM particle mass in this scenario up to 
~ 50 keV. At the same time, for the same mass the primordial velocity distribution of RP sterile neutrino DM is 
colder than in NRP one. This f m ax is as much as the order of magnitude bigger than (|21[) (c.f. [68]). This brings down 
by a factor ~ 2 the analog of the mass bound (|34|) . Analyzing available spectra for a range of lepton asymmetries, 
we see that models with m RP > 1 keV are allowed. Thus, there is a large open "window" of allowed DM masses (c.f. 
Fig|4]). However, as the dependence of the velocity spectrum on the lepton assymetry is not monotonic, to obtain the 
exact shape of the lower bound on the mass at given mixing angle more work is needed. Nevertheless, our results 
show that the sterile neutrinos, produced in the presence of lepton asymmetry, are viable DM candidates, allowed by 
all current bounds. 

Finally, we would like to comment on the mechanism of production of st erile neutrinos from decay of massive scalar 
field, for example the inflaton [6!| (for other mo dels see I 57l . [70l . l7ll . Il55l |). The primordial phase-space distribution 
function for this case was computed e.g. in (69l. [70l . [l55lll56l ]. Maximal value of phase-s pace density for this distribution 
is that of degenerate Fermi gas. Notice that the distribution functions in (69l. l70l . fl55| f(p) is formally unbounded for 
small momenta: f(p) ~ p^ 1 ^ 2 . From this one can easily find that the fraction of particles, having maximal phase-space 
density, is ~ 10~ 8 . As only this small fraction of all particles has maximal phase-space density, we expect the mass 
bound in this case to be stronger than (f3"2"|) . The detailed analysis will be presented elsewhere. 

After this work has been completed, we received a d raft of the paper |l57j ], where similar issues have been considered. 
Our results are consistent with those of discussed in |l57j ] wherever they overlap. 



Most of these studies use the quantity Q(r) = p(r)/cr 3 (r) as a PSD estimator 
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APPENDIX A: INFLUENCE OF ASPHERICAL SHAPES OF DM HALOS 



We analyze the change of the bound eq. |T]) due to the deviation of a DM halo from a spherical shape. Such 
asphericity affects both the spatial volume V and the escape velocity Vx,. We consider the dSph as homogeneous 
ellipsoid with semi-axes a, b and c and assume the ellipticity of its 2D projection 14 e < 0.5. Because we observe only 
2D projection of such an ellipsoid, there are two possibilities: 

Prolate dSph: c > b ~ a. We see the axes b and c, related to the "averaged" radius R via b = R(l — e) 1 / 2 , c = 
R(l — e) -1 / 2 . The spatial volume V is therefore 

V = -nabc « ^7ri? 3 (l - e) l/2 « ^7ri? 3 (l - 0.5c). (Al) 



The gravitational potential for e < 0.5 is dominated by monopole and quadrupole components, 

0«0(°)+0( 2 ). 

The maximal value of the potential occurs near the end of the minor semi-axis: 



5 

„3 



^^oo I spherical 

which gives us the correction for m DEG of smaller than 1% (for e = 0.5). 

Oblate dSph: c ~ b > a. We observe the axes a and c, therefore the spatial volume V changes by (1 
1 + 0.5e. The maximum of the gravitational potential is then given by 



(A2) 



I _ v 1c G N M G N D ZZ 

\4>maA = — = — (A3) 

where D zz = 2M ( C ~ a ) _ the quadrupole moment of the system |l58j . For e< 1 we then obtain 



Vv oo\prolate ~ ^ + 0.05e, (A4) 



-1/2 



i _ G N M G N D XX G N M . 
W ma x\ ~ — — + 2fl3 » — — (l + 0.1e). (A5) 

where D xx is given by the same expression, as D zz above. The maximal phase-space volume changes in the 
oblate case by w 1 + 0.65e, so the correction for m DEG will constitute about 8% for e ~ 0.5. 

Thus, the departure from spherical symmetry for DM halos of dSphs changes the limit on m DEG by less than < 
10% for the case of axis ratio 1:2. This uncertainty is below several others, therefore, we will consider dSphs to be 
spherical in what follows. 



http : //sec . bitp . kiev . ua 

Throughout this paper, we define the ellipticity e in a way similar to that in |88j (see also |83j). i.e. e = 
semi-major and semi-minor axis, respectively. Thus, the case of e = 0.5 corresponds to axis ratio 1:2. 



1 — 6/o, where a and b are the 
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APPENDIX B: ENTROPY FOR DIFFERENT DISTRIBUTIONS 

In this Appendix we will calculate the entropy for several phase-space distributions, including those of ([4]), (|20[) . 
(JSJ), and explore its relation with the quantity Q, defined in 0. 
The entropy of an ideal Fermi gas is given by the expression jl59l | 



S = - I d 3 pd 3 - 



/(r„) log I ^WM^) + (' - f { r, P )) log (l - 



(Bl) 



,(2^ ^■»'>J"°y g 

If the distribution function f(r,p) <C ^n) 3 ' we obtain the expression for the entropy of a non-degenerate ideal gas: 

f(2xh) 3 f(r, P ) 



S = - J d 3 pd 3 rf(r,p) 



log 



<J 



(B2) 



1. Ideal Boltzmann gas 

We start with the case of ideal Boltzmann gas: 

f(r,p) = foe-^r . (B3) 



Substituting it into Eq. (|B2[) . we arrive to the well-known expression (c.f. e.g. [159 . §42]): 



S 5 , fgV (mT) 3 ' 2 . 

'tog ^ LL ■ ( B4 ) 



iV 2 6 V N (2tt^ 2 )3/ 2 

where is the volume of the system, N is a number of particles. Expressing S/N as a function of /? and (v 2 ), we 
finally obtain relation between the entropy and Qin the form (|10j) 



- = logC B - log C B = 57^7^ « 3 x 0.1489 . . . (B5) 

1. Isothermal phase-space density distrubution 

Next, we consider the case when the PSD distribution can be approximated by (pseudo)-isothermal sphere (c.f. (JSJ)): 

, . 9a 2 p 2 , N 

fiao{r,p) = ——— 2 2 , 3/2f 2 , 2 , e ^ T ^ 1 . B6) 

The number of particles in such a system, as well as the total entropy, diverges for large r, however the entropy per 
particles grows logarithmically at large r and therefore the exact value of cut-off is not important. 
Truncating the expression for the entropy at some r max and taking r max 3> r c , we obtain 

|=-log^+logC iso , C». = g SP (1 / 3 % * g x 0-0604... (B7) 
TV to 4 V3(27r) 3 / 2 

3. Entropy for Fermi-Dirac and NRP distributions 

Next, we analyze the case of primordial momentum distribution, which has the form of (rescaled) relativistic 
Fermi-Dirac. 

/(P)= (M)VM^ + 1 ' e(p) = P - (B8) 

For now we keep both F and T to be arbitrary. The distribution in the form (|B8[) accounts for both Q and (|20[) 
cases. The entropy of TV particles with distribution (|B8[) is given by the expression (|B1[) . which reduces to 

qVT 3 

S = — — »I(F), (B9) 
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where function 1(F) is given by 



/>oo 

1(F) = dzz 2 
Jo 



+ 1 \ v e z + l/ \ e z + 1 / \ e* + l. 
The integral (|B10|) can be computed numerically. At F -C 1 the expression (|B10|) can be approximated by 



(BIO) 



3 f 00 dzz 2 

I (F) ^ -((3) (F-F log F) + F __log(e z + l). (Bll) 



The specific entropy S/N equals to 

S_ _ gm i I(F) / C(3) \ 3/2 (v 2 ) 3/2 _ 9 1(F) f C(3) \ 3/2 _m 4 



A 2vr 2 ^ 3 Vl5C(5)y p 2n 2 K^C(5) J Qh 3 ' (B12) 



Therefore, we see that for the distributions of the form (|B8[) relation between the entropy per particle and Q is not 
given by the simple expression PH|) . 

Up until this moment we kept parameters F and T in (IB8[) independent. However, we are mostly interested in two 
particular cases: (i) F = 1 while T = T FD - arbitrary (distribution ([1])); and (ii) F < 1 having arbitrary value, while 
T being fixed to T v - the temperature of neutrino background, related to the temperature of the CMB background 
today via T Uo = (4/ll) 1 / 3 T CMB , (distribution (J2DJ)). 

We start with the case (i). Expressing p as a function of T FD , we obtain 

Q = , (B13) 



where numerical constant q is given by (c.f. [27j|): 

C 5/2 (3) 



20tt 2 v / I5C 3/2 (5) 



1.96... x 1CT 3 . (B14) 



As a result for the distribution (j4|) and fixed number of particles, the quantity Q is independent on T FD , volume or N. 
The entropy per particle is also independent on both T FD and V and is given by 

4 = s = ^(1)tt^tt -4.20... (B15) 
N v 3C(3) v ' 

Although both quantities S/N and Q are simply constants, we find it convenient to choose them in the form (flU)) : 



A = ^ l0g U 



lo g( ) + 1o S C 'fd ! C FD =.g-q-e s «.gx 0.1311... (B16) 



In case (ii) when F <C 1 we obtain for 5/ AT: 



A 3C(3) F ° 7 3C(3) 

Similarly to (|BT3|> - (|B14j) 



2 m ~'l-lo g F) + *. (BIT) 



<?qF. (bis) 



Combining (|B17|) — (|B18|) we can write 

S , /Q^ 3 \ . _ /, 2i 



- = -log(^ r J +logC NRP , C NRP = 5 qexp(l + — ) « fl x 0.137... (B19) 
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dSph 


Qf 






?71fd,hd 


?TInrp,hd 




jfg) ( km 
pc^ V sec 


r 


[keV] 


[keV] 


Sextans 


5.681™* 


■10 


-6 


n O94+0-098 

U.OZ*_g q 57 


1-U4_ 24 


Fornax 


8.86+^7 


•10 


-6 


362+ 091 

U.ODZ_ 051 


1 20+ 042 

i -^ u -0.25 


Leo I 


1 rr+3.08 
r-OJ_o.85 


10" 


5 


416+ 0131 

U.110_g q 75 


1 4^+0-64 


UrsaMinor 


-1 70+2. 70 
1- '8-0.90 


10" 


5 




, r 9 +0.55 


Bootes 




10" 


5 


520+ ' 064 


1 O5+0.33 
1-»J_0.30 


Draco 


n 9-1+1.27 
°' z± -0.82 


10" 


5 


4OO+0-044 


1 85+ - 22 
J - oo -0.17 


Carina 


2.601?; 42 


10" 


5 


n 474+ - 111 

^-0.072 


-1 7 9+0.56 
L - ' -0.34 


Sculptor 


5.93^; 3 3 ° 


10" 


5 




9 9 7 +0.45 
' -0.35 


Leo II 


6.39l°;?° 


10" 


5 


5QQ+0-115 


9 O9+0-62 
^■•J^-0.40 


Canes Venatici I 


6.16it:S 


10" 


6 


n oon+0.045 


i - u ' -0.13 


Ursa Major I 




10" 


5 


44n +0 ' 077 




Hercules 




10" 


5 




1 74+0-67 
x - '^-0.43 


Leo T 


6.26±|;I| 


10" 


5 


0.590lg;SS 


9 oi +0.85 
^"Jl_0.48 


Ursa Major II" 


-, in+1.61 

i - u -o.5S 


10" 


4 


0.6851S;1S 


9 oi +0.96 
^•81-0.55 


Leo IV 


o or+40.91 
2 27 


■10 


-4 


0.8981°;^ 


4.04l5; 5 27 


Canes Venatici II 


s.git^ 1 


•10 


-4 


i rn +0 - 30 

J - uo — 0.18 


4.88li;? 


Coma Berenices 


5.4611;?; 


10" 


4 


1 ni+ - 28 


4.751J; 82 



TABLE III: The mass bounds, based on the evolution of the average PSD Q [26l. |27|. The bound is provided for illustration 
purposes only (see Section [C] for discussion). 



APPENDIX C: MASS BOUNDS FROM THE EVOLUTION OF THE AVERAGE PSD 

For illustration purposes we provide in Table HTT1 the average PSD estimator Q for all the dSphs, considered in this 
work, as well as the lower mass bounds, based on the inequality © for Q during the evolution [2|| [27| (for detailed 
discussion see Section [TTJ) . The value oi Qf, shown in the second column of the Table HTT1 is calculated from the data 
in the columns (3-4) of the Table [Til using formula (|26[) (with rj = 1) and Qi is defined via ([5]) for the momentum 
distributions ^ and (|20[) (for the bounds 7ti F d,hd and rn NRP HD correspondingly). The results for Leo IV are quoted 
in ||5BJ) (Section E| . 
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